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ABSTRACT 


An airborne tactical system composed of three 
distributed microcomputers is described. An 
extrapolation technique using the method of orthogonal 
polynomials is presented to solve the ballistics 


problem. 


In order to test the performance of distributed 
microcomputers for real time military appiications, a 
Simplified airborne tactical system, utilizing two 


miccocomputers working in parallel, is implemented. 


Themiem=varnon £Lor chis tresearch ts the reduction 
Or cost that would result from the use of 


microcomputers in such systems. 
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Airborne tactical systems on board attack aircraft are 


designed to rferform the following functions: 


1. To navigate to a given target. 


poet Oo COnNtionwmene: auto-pilot. 


3. To display information to the crew members. 


YW. To solve the bailistics equations that will 


compute the down range travel of the weapon. 


De lomeecGerapotate rrom the available data, the 


release point of the weapon. 


Current tactical systems employ a minicomputer as their 


main computing device. 


An alternate approach to solve this problem is by the 
use of several microcomputers processing in paraliel in 
place oF the minicomputer. Although the most recent 
microcomputers are much faster than their predecessors, they 
still lag behind the minicomputers in their arithmetic 
processing capability. That is why a distributed system is 
needed. In such a system, the processors are put to work 
with as much parallelism as possible and they communicate 
With each other by means of an interrupt mechanisn. The 
result is that the single conputing cycle is broken up into 
Smaller cycles some of which are executed in parallel, 


reducing the response time. 





In this system, three microprocessors constitute an 


airborne tactical system: 


1.The Executive computer. 


2.The Navigation computer. 


3.The Ballistics computer. 


Because of the availability cf only two microcomputers, 
this thesis implements a Simplified version of an airborne 
tactical system using two Intellec/8-iNod 80 microcomputers 
working in parallel to simulate the attack function of tne 


system. 


The system computes from a given series of aircraft 
positions with respect to the target the release point of 
the weapon so as to hit the target. All of this happens ona 


real time basis, controlled by an external clock. 


Fig 1 shows a photograph of part o£ the systen. 
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To provide the reader with an understanding of how a 
tactical system performs, Chapter II describes a_esysten 


uSing three microcomputers. 


Chapter III consists of an analysis of the ballistics 
problem which will be the critical part of the implemented 


systen. 


In Chapter iV, the work done on the I.B.M. System /360 
With FORTRAN is described. The objectives of this phase of 
the Bee ject were Mainly experimentation with system 


Organization and design and testing of numerical methods. 


Chapter V haS a brief report about the necessary 
hardware adjustments and setup needed for the 


multiprocessing scheme. 


Finally, in Chapter VI, the PL/M implementation is 


described, together with four test cases. 


The data used were partially obtained from real runs 


performed at the China Lake target range. 


Appendix A contains the FORTRAN program that wasS written 
to perform the converSion from the encoded data of Lange, 
elevation and azimuth of the aircraft, as tracked by a 
Nike-Zeus radar, into decimal representation of altitude and 
horizontal distance to target. The program also computes 
speed and dive angle, to be used to test the extrapolation 
routine. This same program, converts these values into the 
three byte floating point number format for the PL/M 


vrogranh. 


Appendix B shows a program which tests the extrapolation 


rouciane. 
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The system 1S composed of three computers, sensors, 
display units and auxiliary devices, each of which is 
described below. The reader may find a more complete 


description cf the sensor instruments in [2}. 


A. THE NAVIGATION COMPUTER 


This computer executes a program that periodically polls 
sensors in order to compute the change in che aircraft's 
position during the last tine increment.There are four sets 
of sensors to gather the necessary analogue data which is 
then converted to digitai form and presented to the 


computer.These sensor instruments are: 


(fimemiiertial Navigation System (INS). 


Z-ihe Doppler Radar (DR). 


Setnewnir Data Computer (ADC). 


4&.The Radar Set (RS). 


The INS is Zi gyro~stabilized platform with 
accelerometers mounted so that accelerations of the vehicle 
ace measured in the north-south, east-west and vertical 
directiCOnsvweesumoverates on the principle that every time the 
aircraft changes speed or direction it experiences an 


Le 


acceleration. With the continuous measurements of 





aircraft's accelerations the INS determines the velocity by 


means of analcgue integration. 


A Doppler Radar radiates a pattern of beams to the 
surface of the earth and receives the reflection of this 
energy back. The difference between the frequency of the 
transmitted signal and the frequency of the received signal 
caused by the relative motion of the source of the signal 
With respect to the reflecting surface is known as _ the 
Doppler Effect. By measuring the frequency shift along each 
beam, the DR computes the velocity and distance travelled by 


ne alrcraft. 


The Air Data Computer consists of the aerodynamics and 
thermodynamics sensors: angle of attack vane, static 
pressure source, Pitot tube and temperature probe. [It 
computes the airspeed, the free-stream outside-air 


temperature and the Mach nunber. 


In order to use the ADC for dead-reckoning, the attitude 
and heading cf the aircraft must be supplied from other 
sources. This information together with the information of 
wind's velocity and direction, allows the determination of 
the aircraft's position by extrapolating from a previously 


known position. 


The Radar Set determines the distance of the target by 
measuring the time elapsed from the transmission of a 
electromagnetic wave to its reception after being reflected 
by the target. The bearing is determined by the antenna's 


Grrece lone 


The system has a magnetic compass that indicates the 


aircraft's heading in relation to the magnetic north. 


An altimeter supplies the systen with altitude 
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Pmtorna tion. 


There are four basic modes of operation of this’ systen: 
By using both the INS and the DR, the system operates in its 
most accurate mode. If one of these primary sources is 
malrunctioning, the system operates without its information, 
resulting in other modes of operation. In the case of 
failure of both INS and DR, the system operates in its least 
accurate mode using the ADC with the estimated wind vector 
and true airspeed of the aircraft to compute latitude and 


longitude. 


No matter in which mode the system is operating, the 
information of aircraft's position is smoothed by the 
Navigation computer. The four most recently computed roints 
acre approximated by a quadratic polynomial uSing a 
least~squares fit. This smoothing process has a damping 
effect that reduces spurious oscillations of the data 
points. The mathematical description of this technique is 


heund in Chapter III. 


Using the smoothed data values, the increments in 
distance travelled in the last time increment are computed 
and the position coordinates displayed. AS soon as the new 
values are computed, the Navigation computer interrupts the 
Executive and transfers this information to the Executive 
computer. The Navigation computer has the ability to 
function in a stand-alone mode to display aircraft location 


regardless of the status of the rest of the system. 


Another important feature of the least squares curve 
fitting technigue is the ease with which sensor malfunctions 
May be spotted. The crew is warned of the ill behavior of 
one of the sensors, and the system's wode of operation may 


be changed. 
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fens EPADLISTICS COMPUTER 


This computer receives from the Executive computer: 
aircraft speed, dive angle and altitude. The type of weapon 
and target's altitude are initially received from the 
Executive computer. With all these data, the computer 
numerically solves the differential equations to compute the 
down range travel of the weapon. There will be no guidance 
of the weapon after release, although the weapon may he 
rocket assisted or have a retarded mode of fail. When the 
Ballistics computer finishes its calculations, it interrupts 


the Executive computer to transfer the results. 
See fLHE EXECUTIVE COMPUTER 


The Executive computer controls the function of the 
whole systen. ifesenae=eo the Ballistics computer, at the 
beginning of the run, the target's altitude ana the weapon 
selection. The Executive computer will be interrupted by 
MIomENavrgatz on | computer which periodically reports the 
eEecrake Suu position. This information is sent to _ the 
Ballistics computer after it is converted into the proper 


fOr Ms 


The number of data points that the Executive computer 
needs for the extrapolations depends on the degree of the 
Pitnigueueve. if this system, it will be a third degree 
polynomial. In order to guarantee a solution to the 
extrapolation process, the Executive computer has to be 


Supplied with at least four distinct data points. 


When the Executive computer has enough data points 


representing horizontal distance to target and respective 


15 
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down range travel together with the corresponding time, it 
extrapolates these two functions o£ time and finds their 
intersection. At the intersection of these two curves, the 
distance to tne target wili be the same as the down range 
travel. Thus, this will be the release point which is 


Sought. 


fies cemeputer will divide its time between data 
fiorpiulation, data transmission, displaying of infcrmation, 
computation of release point and delivering steering 
commands to the auto pilot. If so desired, when the release 
point is reached, the conputer will automatically trigger 


the release mechanism. 


An external real tine clock provides the Executive 
computer with exact time, so that time dependent variables 


may be recorded together with real time values. 


ieee tis CCHMUNTCATION SCHSYE 


Because of the inherent hierarchy between the computers 
involved, the master-slave type of multiprecessing is the 
most suitable and simple form to implement. The Ballistics 
and the Navigation computer act like peripheral devices of 
the Executive computer, resulting in a one way interrupting 
scheme. The only computer which has to be interrupted is the 
Executive computer. The Navigation and the 3ailistics 
computers are the dedicated slaves which asynchronously 


interrupt the Executive computer. 





III. THE BALLISTICS PROBLEM 


In the most general form, this problem may be described 


as follows. 


The aircraft's position in space is described at time t 
by a radius vector r(t) in a three dimensional coordinate 
systen. The aircraft's velocity is described by the 
derivative with respect to time of the radius vector, i.e. 


mmenvector dr/dt=r 


These two vectors are shown in Fig 1 
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It is assumed throughout that the aircraft during its 
final approach to the target is restricted to fly in a 
vertical plane that contains the target and the origin, so 


that the problem becomes two dimensional. 
In the presence of strong winds, the ground track of 


aircraft is used to determine the X-Y¥ plane as arg 3 


shows. 
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The target 1s located on a line which is parallel to the 
@eound track Of the aircraft, so that at the estimated 
release altitude the wind will carry the weapon thrcugh the 
distance between the ground track line and the dotted line 


Baballel to the xX-Y¥ plane. 


gn the case that the wand direction coincides with the 
aircraft's heading, the target will be in the X-Y plane and 
mmemwina effect will either reduce or increase the down 


range travel of the weapon. 


With the problem, now reduced to two dimensions, the 


foliowing mathematical formulation applies: 


At the instant the weapon is released from the aircraft, 
the initial position and velocity of the weapon agrees with 
mies pOs2tion and velocity of the aircraft. There may be a 
small added initial velocity when the weapons are released 
by small charges, cr there may be a large added muzzle 
Metecity £er rockets and guns. In any case, the initial 
velocity of the weapon is determined by the aircraft's 


velocity EE ieplus the constant vector V. - 


After released, the weapon becomes a freely falling body 


satisfying Newton's differential equation 


mu=mg-c Jul u, 


where u is the position vector of the weapon ana the 
forces acting on the weapon are the gravity and drag which 
PoP eoreomtmonal to the velocity squared. The term c is a 


characteristic drag term and it 1s weapon dependent. 


In the two dimensional coordinate systen, the equations 
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become u=(x,y). 


" ar 
eae k +y <x 


s: eames oe 
my=-mg-c \/x +y y (1) 


Mech inetlal conditions: 


oa) e— t ia =c ot 
x ( i ra a X ( 9 a ae 


t )=r (t Pies jyesemats 
ao an 9? es . ue 


The bembardier's problem is to determine the time t at 
Wmich to release the weapon in order to hit the target. 


Stated mathematicaily, what needs to be _ solved 


combination initial value-boundary value problen: 


mx=-c\/x +y x Meet yer {t) X(t) =C (t} tv 
x x 


» 





p o.2 Z 2 Steet coe 
jee: y y(ty=r_ (t) y(f)=r (bt) tv 
Y Y y 
With the boundary condition that at time t: 
y(t}=0 x(t)=? 

In order to give a clearer understanding of the problen, 
let it be further simplified by assuming that the diag 
forces on the weapon are negligible. In this case, the 


differential equations become very Simple. Let it also be 


that the release point an The system becomes: 


pM 





nx=0 a= tO x10) =c 9Cey tv 
Xx x XK 
my=—-ng y(0)=r (0) y(0)=rc (O0)+v 
y y y 


The analytic solution to the initial vaiue problem is 
obtained by integrating both equations twice with respect to 


t, resulting: 


X(tj=(c (0)+v )t+r (0) 
a4 x »6 


2 ; 
V(t} -sge 72+ (te (0) tv )tt+r (0) 
y " y 


en Py : ao ; : * 
Mest 1b es rdesaeca tO Lind at™ What time t does the 


; , BO acs , : +, , eee 
mempon hit the greund, 1.e. at wnat tite t 15 the condition 


mace) = 0 met. 


This question is answered by solving tne quadratic 


equation for the positive root: 
2 é 
-gt /2+(r (O)tv )ttr (0) =0 
x y Y 


giv wage 


es i 2 
m—(e (O)+v )7gt \/({e (0) tv) +t29v_ /9 
y y Y y Y 


Therefore, in order that the weapon hit the target at 


time t: 
T=x(t)=(c£ (0) tv )ttr (0), or: 
x » 4 Xx 


i 58D 
f=4E (0 r {O)tv )faot\f(c (O)tv ) t2gr /g}tr (9) (2) 
(E (ev ILE O49) gt VE (0) tv, , x 
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hoasan Cxample, af the aircraft is at level flight with 
constant speed: 


ay a? 
Bare Ge) V +20h /gj+r (0) 
x y y x 


bao Sreiease point (xr (0),h) is then given explicity by: 


» 


[DP 
Beep =t-—(Stv Viv /gt\jv +2gh /g) 
x im y y 


Even in the general case, relation (2) allows the 


determination of the release point (r (0),r (0)) itplicitly. 
X 


Therefore, if the drag effect is neglected, the 
initial-bcundary value problem is relatively easy to solve. 
However, drag cannot be ignored in the real case, requiring 
that the more complicated system of differential equations 
(1) be solved. 


If a dedicated computer is employed to solve the initial 


value problem (1) sufficiently fast, the down range travei 


BE may be determined for any release instant a : 


Therefore the release time t would be that value for 


Mich x(t)=T. 


To solve this problem, all that is needed is an 


extrapolation technigue, so that from successive aircraft 
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horizontal distances to target and corresponding computed 
down range traveis it 1s possible to estimate future values 
for these two functions. At the release point these 


estimates should coincide. 


The extrapolation technique employed is a least-squares 


approximation as follows. 


It is desired to approximate a series of N data points 
SE (x } by a function of the fora 
n 


ees eet (XK) tact w(x) s 


that depends linearly on the parameters CP aoe 


tice tuneutonus U0 (x), are selected a priori and the 
u 


parameters must be determined. 


(Peerderss to approximate the data points f(x ) by the 
n 


function F(x), the parameters must mininize the function 


te 2 
mde send y= ) [ite j= (x) 7 
1 k nz nh 


eremeecriterion function teads to the least~Squares 


method. 


By applying minimization techniques, one gets the 


following system of k linear equations in k unknowns: 
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= * 
Cee (cy Ue ke) +t. +d Lu (x yu (x )= 
hn 1 eR kKns 2 n k on 
N 
=) £(x )u (x) a= 1 ee ok (3) 
net ioe i fh 


It may happen that although this system always has a 
solution it turns out to be unreliable if the matrix A of 


the coefficients is ill-conditioned, i.e. if its condition 


number, cond(A)=|{Af{ JA | is too big. 


Next it will be shown what the effects of using 
orthogonal polynomials are and how to generate these 


polynomials. 


The cond(A) depends entirely on the choice of the 


mimetions wu (xX). One way to garantee a smali cond(A), is to 
at 


use orthogonal functions, to get: 


du. (x yu (x )=0 1#j 
sell 0) ey Gig eg 


~<— 


By such a choice of orthogonal functions the system (5 


reduces to: 


x «(NM ; N 
d. you (x ju (x)= L £(x Ju, (x } =p ecewear 

Mor! 2 ft nn Ne” fee 

which, besides rendering reliable results, is much simpler 


to solve, reducing the computation involved. 


The definition of the scalar product of two functions g 


and h that will be used is: 
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N. 
<g,ho= L g(x )h(x ) 
mel n n 


The functions are said to be orthogonal to each other if 
<g,h>=0. 


A sequence of polynomials (2, G9] 1s said to be 
Smcnogonal if all ee) are orthogonal to each other and 
each eo) is a polynomial of exact degree i. One property 
of such a Sequence is that any polynomial p(x) of degree <k 


can be written as 


= 2 Cee) fe. 5 51d 2 tt 
Pp (x) oh mo : 5 as 7 k ee (4) 


The measure of difrerence between £(x) and p(x) that 


needs to be tinimized can be written as 


i 2 
<£ (x) -p (x) -£ (x) p(x) >= DD CE(x )-p(x ) 1] 


And, 


)O (0G) a he) 
sa ao) 


re - ~eeen ne -~d P -~...-d P > 
SEX) aco Sn Se (x) ; il : 18 


» 


The coefficients that produce the best fit have to 


satisfy the system of linear equations: 


2 Ps * 
Gomera <P (P >t...¢d <P ,P >=<£,P > L=0 owes 
Ome 07) 3: 1 74 G2 k K .L 1 





PemOrELHOGOnality, this reduces to: 


Gute eee =0,...,k 
aE 1 


the ccefficients will be given -by: 
ed 

dmea><t,P >/S. i=) eee 3 
1 i 1 


tow 1it will be shown aow the sequence of orthogonal 


polynomials may be generated fron the N data points. 


These polynomials satisfy the three term tLrecurrence 


relation: 


om Cpe temoces YF (x)= C P. (Xx) V0 7 fee 
1+1 2: hg ee 1 i-~1 


Beceem (x) ,P (x) >/S. Tas Kod 
at 1 a i 


C =arbitrary 1=0 
1 


mony (A oS}  i>0 
1 11 1-1 1-1 


Where A =a “/2@. ali 1 
mec t. 2 


S =<P.,P >#0 = Greener ko and a. is the ieading 
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coefficient cf the polynomial P 
i 


1s decided to get all orthogonal 


HOG Sitplicity it 
pelynomials with leading coefficient equal to 1, so that: 


aly. 


Setting Ee 


N 
S =<P ,F >=) 1=N 
0-0. Oo 1 


Next, from the recurrence relations: 


Pe (X)=(x-B )P (x)=x-B , where 
44 ee a o ! ) i 


WN 


IOP ek ak (xX) O/S = » x 5 
Q pee syle ican ae: 


With ce 2 (X) already obtained, the general step 


is: 


N=4 


N 
te Z 
ee ces 2, [P(X ) J 
3 J 
1S ©f exact degree 4, 5S can be zero only 
J 


Because P (x) 
J 
distinct.) oc 


if no more than j of the points x ince are 
there are more than j Gistinct points among the data points, 


it is possible to calculate: 


Ba oxXE) (X),P (xX) >/5 . 
J 2 2 ek 


eo 





7 
3) a 


The next orthogonal polynomial will be: 


EC (ee (x}-C P (Xx) 
j+1 2) ae J 


As an example, the method will be applied to approximate 


the second degree polynomial 


2 
fee) —10—2xX+x 710 
Moeete pPCantsS Given by: 


x =10+ (n-1) 75 ft =f(x ) Resihgteters-7 0 
n n n 


¢ 
2 

ae ~p(x )j] , needs to be minimized. 

La | n 


Thus, 


6 
P (xsl, See nes 
0 { 


G 
and B= >» [10+ (n-1) /5 ]/6=10.5 


rz 4 
& 2 
as i 
P (x) =x-10.5 a pi @-1) 75-0.5 3} =0.7 


nei 


G 
2 
B= » { 10+ (n-1) /5 Jf (n= 1) /5-0.5] /0.7=10.5 


n= 


€ =s /S =0. 1166667 
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2 
P(x) = (x- 10-5) -~0. 1166667 


* Ox * 
Next the coefficients aS and A for the least-squares 


approximation: 


x x 2 x 
x)=d P (x) +d P (x)td P (x are calculated: 
ce eo Fe Ne a 


te G 
Me 60 102666667 
O mn n 


b 

d = Lf Patx 7 0. 7=0..1 
nm niin 

me Ls 

a = ).£ P (x )/0.05973332=0. 999939 
nt Tf 2 R 

Finally: 


* 
pee — 0 - 0 G0 0060 710. 1 (x- 10. 5)=t 


2 
Ome (kak. >) —0.1156667 } 
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IV. FORTRAN SIMULATION 


et cee ee ee = ee Se ee ee eS eS 


It Was realized that it would be advantageous to produce a 
FORTRAN Simulation using the IBM System/360, because: 
1. The computing environment has more power and the 
acithmetic is easier to handle than in PL/SH. 
2. The system's performance can be checked with 
respect to timing constraints. 
3. The alternative computer program allows accuracy 


comparisons to pe made with the PL/M implementation. 


After the FORTRAN simulation had been tested and fully 
exercised, the PL/M program was then produced, taking 
advantage of the experience obtained during the development 


and use of the FORTRAN progran. 


In this Simulation, the Navigation and the Ballistics 
computers are represented by subroutines and there has been 
no attempt to simulate the transfer of data mechanisn 
between computers, because it was realized that the aznount 


of effort reguired would not justify it. 


The modules that comprise the simulation and their 
PUnNGtPOnsS are: 
HAIN Governs the togic of the simulation and takes care 


of the real time clock. 


INIT Performs the initialization for each particular 
Ut. 
ICHECK Checks if there is a solution to the problen. If 


there is a solution, it then verifies 1£ there is 
enough time to compute the solution. If tke 
aircraft is too low and could be hit by the 


weapon's fragments, the pilot is warned. 





DSPLAY 
POS 


BAL 


SETDAT 


DECODE 


TRAJ 


RUNCE 


DERIV 
INTER 


Cou P 


Cher tT 


Displays the resuits and messages. 

Galeulates the aircraft position and updates its 
airspeed. 

Interfaces the MAIN module with the Ballistics 
modules. 

Does the initialization process for the Ballistics 
computations. 

Sets parameters for the Ballistics computation, 
according to the selected weapon. 

Caiculates the downrange travel and time of flight 
of the selected weapon at the current aircraft 
position and speed. 

Auxiliary module of TRAJ. 

Auxiliary module of TRAJ. 

Commands the intersection search pattern for the 
horizontal distance to target and the downrange 
travel curves. 

DoeS a stepwise comparison between the horizontal 
distance to target and the downrange travel 
curves, searching for their intersection. 
Extrapolates values for the horizontal distance to 
target and the downrange travei curves for a given 


time. 
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PROGRAM VARIASLSS: 


VARIABLE DEFINITION 

A Runge Kutta parameter. 

AA 0.5/Runge Kutta parameter. 

ALT Weapon release aititude. 

AN1,AN2 Runge Kutta variables (X direction) 

AP1,AP2 Runge Kutta variables (Y directicn) 

eC Meer rExeon drag curve coefficients. 

or GrvGeaeretch factor. 

eCrokit 1 SiiicusoteetGn Lactor for the first stage. 

CrORN2 CKDG stretch factor for the second stage. 

CKDG Bonb coefficient times drag coefficient. 

CM Mach number, 

eT Netnix on Mach Gurve coefficients. 

DEG Weapon release angle. 

DEL ArcTan (ejection velocity/weapon velocity} 

DKG Simm) 1n CKDG 

DKG1 Spit un CKDG Lor the first stage. 

DKG2 Shift in CKDG for the second stage. 

DM Shift in Mach number 

DMAX Largest integration step size. 

DM1 Shift in first stage Mach nunber 

DM2 Shift in second stage Mach nuoker. 

DS Integration step size. 

DL Piberdiotibedratlone step SiZ¢. 

FN Polceete tne thtust £or the first stage. 

FRACT Parameter for Runge Kutta integration 

G Acceleration due to gravity. 

HH Lordi drpeg function: 

LBOTH Key determining if two drag curves are 
reguired by the weapon. 

IDNO Tdentification number of each weapon. 

Sea re Reference for Mach curve coefficients. 

IREG Mach region index. 

MSTG Weapon stage index. 


34 





RAD 


YO 

ba 
NAS 
HV 

VV 
HVA 
VVA 
XAT 
hE yBo 
Peis 
HTINE 
DTEAL 
DTE 1 


DTE2 


DE 3 


Conversion factor -for radians. 

Air density. 

Slope coefficient for retarded snakeye 
fins deployment time. 

Elapsed time of flight from weapon 

release. 

Thrust/mass of rocket. 

Adjusted weapon release angle. 

Alrcraft speed in feetSeconds. 

Velocity of weapon. 

Component of V in the X direction. 

VX at release point. 

VX at start of present integration step 

Component OF V an the Y direction. 

VY at weapon release. 

V¥ at start of present integration step 
Muzzle velocity of guns. 

weapon ground range from release. 
Weavon altitude above sea level. 

Y at start of present integration step. 

Target altitude above sea level. 

Speed of aricraft in knots. 

Horizontal component of VKTS. 

Vertical component of VKTS. 

Previous value of HV. 

Previous value of VV. 

Horizontal distance to target. 

Real time. 

Saved value of TIME for DSPLAY. 

Previous value of TIME. 

Time spent by the Bailistics Computer. 

Time spent by the Executive Computer to 
transfer data from the Navigation Computer. 

Time spent by the Executive Computer to 
ransfer data from the Ballistics Computer. 


Tine Spent by the Executive computer to 
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compute the release point. 


DTEY Time spent by the Executive Computer to 
execute check number two. 

Dir 5 Time spent by the Executive Computer to 
transfer data from the Navigation Computer 
and execute check number four. 

DTE6 Time spent by the Executive Computer to 
save and restore values when calling the 
Navigation Computer and execute check 
number three. 

£C Time necessary to compute for the first time 
the release point. 

WARP (2) Time at release point. 

XARP (2) Horizontal distance to target at release 
point. 

AARP Aircraft alititude at release point. 

THV Time it takes for the aircraft to start the 
evasion. 

ATOL Limit angle between aircraft path and 
line of Sight of the target. 

RLIU Radius of a hemisphere defining a danger 
zone around the target. 

ALIN Peeeengediticuae,:Or the aircraft. 


The units of the variables are: 


Distances: Nautical niles. 
Altitudes: Feet. 

Angles: Degrees. 

Time: Seconds. 

Speeds: Knots. 

Piva omOmeOrsAL, thi, DIEZ, DTES, DTEL, DIES, DIE6, 
va ATOL, RLIM and ALIN were assumed Teens 
Simulaticn.Mcst of them are weapon dependent. 

This pxrcegram used a least-squares fit of the third 
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degree, by means of orthogonal polynomials [3]. After the 
fourth aircraft position and corresponding downrange travels 
are computed, the Executive Computer proceeds to Compute the 
first release point data. From this time on, for every new 
position and corresponding downrange travel received, the 
Executive Computer determines an updated release point using 
the four most recent aircraft poSitions with their downrange 
travels. This process goes on until there iS no more tine to 


compute the next release point. 


The way the release point is found is by a stepwise 
comparison between the extrapolated values of the horizontal 
distance to target curve with the downrange travel curve. At 
the release point, these two curves intersect. After the 
curves have been extrapolated from the first four points, 
the compariscns between the extrapolated values of the two 
functions are made from point number five on. When the 
intersection point is passed, the direction of the search is 
reversed and the time increment is reduced. The search 
continues until the difference between the curves is smaller 


than a fixed limiting valve. 


Hag G& showS graphically how the search for. the 


Patersection is nade. 


3) 
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The ICHECK module has 4 entry points and its function is 
momeneck for exceptions. First, with the initial data, it 
checks if there is a soiution to the problem. There are two 
Situations that would not have a solution: 

1. The aircraft dive angle is such that the impact of 
the weapcn would always be short. This happens if, for a 
given weapon, the angle between the aircraft path and the 
line of sight of the target is less than a certain 
threshold. 

Zoe iewedirerare iS not diving and the present 


distance to target is smaller than the downrange travel. 


There 1s another situation in which the solution. should 
not be computed, because at the release point the aircraft 
would be belcw its limiting altitude. In this case the 
SByotem would also display “NO SOLUTION TO THESE DATA". 


Because the computation of the release point revuires 
approximately one second or real tine, it is necessary firs 
to check if the extrapolation can be carried out in time to 
release the weapon. If this can not be done, the system will 
iesoray NOT EMOUGH TIME TO COMPUTS"., 


Every time a release point iS computed, the ICHE 


OQ 
(bt ON 


module checks if at the moment of weapon impact the aircraft 
is inside a danger zone which is defined by a hemisphere 
around the target whose radius depends mainly on tne weapon 
ects tn tris ts the case,the system will display: “ATACRAY 
i>) EN DANGER ZONE". 


After each release point calculation the ICHSCK module 


verifies if there is time to compute a new release point for 


the next position. 
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Ve. THE MICROCOMPUTERS HARDWARE IMPLEMENTATION 


=a oo — ee ee ee ee ae Se ee ee mm tc me ee ee es 8 er ee 


The hardware implementation of the system required that 
the two computers have a second input-output board 


installed, in order to provide all the connections needed. 


To make the system Simple, it was devised so as to use 
the fewest number of input-output ports possible. One of° 
the key facts that allowed reduction of the communications 
paths among the computers was the predictability of message 
lengths throughout the complete cycle of system operation. 
The result is that there is no necessity of a status port 


during transmissions of data. 


The Intellec microcomputer, aS it is furnished by the 
Manufacturer, does not have available externally the eight 
dines that constitute the interrupt port, nor does it have 
an external interrupt request line. It was necessary to set 
up a board to be installed on the Executive computer, to 
bring out these lines, so that the necessary connections 
could be made. Fortunately, this turned out to be the only 


change necessary to use the interrupt feature. 


Another addition to the Executive computer was the 
external ciock with an eight bit binary counter that was 


connected to one of its input ports. 


Pag 5 shows the system's components. In this figure, 
the numbers correspond to the eight bit input-output ports 
of the computers. PTR stands for paper tape reader and DU 
stands £oOr disk unit. In the Executive computer, the 


Mime Ewoe ee anSteucti1On port (11) and the Interrupt Request 
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line (Int) are shown. The Executive computer will be 
interrupted by means of the Restart 3 instruction, which has 
the bit pattern DF. 
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Figure 2 = MICROCOHUPUTERS HARDWARE IMPLEMENTATION 





The system's software is composed of two separate 


programs, one for each computer. 


The program that runs on the Ballistics computer has two 
procedures to take care of the communication with the 
Executive computer. The first one is to receive the weapon 
identification number and target altitude at the beginning 
of the run. The other is to receive speed, dive angle and 
altitude of the aircraft and, subsequently, compute the down 
range travel. The module that computes the down range travel 
was the subject of a previous thesis [1] and is used in this 
system with only minor adaptations. When the computation of 
the down range travel is finished, the Ballistics computer 
interrupts the main computer in order to send back the 
results. If there is already data for the next computation 
of down range travel, this same procedure will receive the 
new data and proceed to compute. Otherwise the Ballistics 
computer will exit from this procedure and return to the 
main progran. In this fashion the system guarantees that 
the Ballistics computer will be used at itS Maximum. This 
policy had to he adopted in view of the fact that the 
computation done in the Ballistics computer is the most tine 


consuming in the systen. 


The system works in such a way that the Executive 
computer will only send data directly to the Ballistics 
computer if it is idle, in which case the Ballistics 
computer will be waiting for data, no interruption needed. 
If the Executive computer has data to send but the other 


computer is busy, a condition that the Executive computet 





knows by testing a flag, the data is sent to a buffer. 


The interrupt service procedure in the Executive 
computer is the one tnat upon completion of receiving tne 
computed value from the Ballistics computer, checks iff the 
buffer has data yet to be processed by the Ballistics 
Semputer. If this is the case, it will take care of 


transmitting it to the other computer. 


The program that runs in the Executive computer is 
composed of the following moduies, with a brief description 
Seeeneir functions: 

Hon: Peicsmibne program to the I/O facilities 
of the CP/M system if no returned 
value iS expected. 

BIONKZ: ¢ fhe same as above, except that a returned 


value is expected. 


Read: Reads characters from the CRT into a 
buffer. 
Penne Diepiavsetue butter starting at the 


specified address. 


Printchar: Displays the specified ASCII character. 


Crit: Sends carriage~return-line-feed characters. 

Getc: Points tes the current character in the 
buffer. 

Readc: Reads one character from the buffer. 

Readp: ReadS a pair of characters from the buffer 


and returns the corresponding decinal 


value. 


All of the above modules are needed to use the CRT as an 
I/O device by the CP/M system [4}. 
Conv: To return tne decimal value of the number 
character. 
Trans: Returns the ASCII character. 


ClOCK: foster eter values Of the external clock. 
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Realt: Displays real time. 

Relpt: Computes and displays the release point. 

Sendbuff: To send data to the buffer. 

Dsplv: Displays the variable value. 

Display: Displays aircraft's speed, dive angle and 
Qetleliae. 

Pow in: Computes the coefficients and constants 
used by the fitting plynorial. 

Peta Floats the byte values of time. 

Search: Finds the intersection of the curves being 
extrapolated. 

Round: Chooses the closest time to the 
Pirersecet lon point, in 
decimals of second. 

Hee cic Given the time,extrapolates the value of 
the function. 

Negat: Tests if the variable is negative. 

Rane. Converts the mantissa of the number to 
decimal value. 

Poten: Computes the decinal value of the exponent. 

Integ: Displays the integer part of the values. 

Decin: Displays the decimal point and one decinal 
place fraction of the variable. 

Index: Finds tke location of data in NAVTAB given 
the time. 

Comp: Auxiliary module of Search. 

Intserv: Interrupt service procedure.’ 
Receives computed values from the 
Bellistics computer. If the 
buffer is full, it sends new data 
to the other computer. 

The Sinulation of the Navigation computer eS 


accomplished by 


the use of a tabie vith valves of speed, 


dive angle, altitude and horizontal distance to target at 


points in time one second apart for the first ten points and 
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one tenth of a second apart for the last thirty points. 


The first part of the table is used to obtain points to 
extrapolate values of horizontal distance to target. From 
the calculated values of down range travel using these same 


points, extrapolations of this variable are made. 


The second part of the table has finer time increments 
in order to allow a more precise determination of aircraft's 


position at release point. 


The complete cycle of operation of the system is as 

follows: 

ir ieccmeuccning 2f the Ballistics computer is 
ready, the Executive computer waits for inputs of run 
identification number and weapon identification number by 
Cal . 

Ze The Executive computer sends the weapon 
identification number to the Ballistics computer. 

3. The Executive computer waits until the clack reads 
zero seconds. 

4G. The Executive computer brings from the NAVTAB 
table to a working area, the data corresponding to the 
present time. 

5. These data are sent to the Ballistics computer if 
memeroetdale OF to a buffer if it is busy. 

6. The Executive computer displays the data sent, on 
the CRT. 

7. The Executive computer waits until two seconds 


after the last point was sent to the Ballistics computer. 
Steps 4, 5 and 6 are executed three more times. 
During the time these events are taking place, the 


Ballistics computer calculates the down range travel and 


interrupts the Executive computer in order to send the 
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results of its computations. 


There will be two instances where the Executive computer 
interrupts will be disabled: during buffer accesses and 
during displaying of data on the CRT. 

8. When the Executive computer already has the down 
range travels for all four points, it computes the release 
point. 

9. With the computed time for the release of the 
weapon, the NAVTAB table is accessed in order to obtain the 
corresponding position of the aircraft. 

10. This data is sent to the Ballistics computer. 

11. The Ballistics computer calculates the down 
range travel and sends it back to the Executive computer. 

io eecitcemeciit 15 displayed on the CRT. 


The system returns to step 1. 


The two second interval between the data points used for 
the extrapolation was adopted because cf the displaying of 
information on the CRT during the time the computers are 
computing the release point. fhis display function takes on 
the order of one second to be _ performed. in ae eeak 
implementation, the displaying would be avoided during the 
release point computation, because it would slow down the 


process. 


This program was exercised with four sets of data. kuns 
one and two used real data that was obtained from a digital 
recording of radar tracking during A7 aircraft practice 
bombing runs at China Lake. Run three used data generated 
to represent level, constant-speed fiight. Run four used 


constant dive angle and speed. 


The sequence of displayed messages cn the CRT, for runs 


one and three, were photographed. These photographs appear 
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in the Section entitled PL/M program output, starting on 


page 98. 


In the fcllowing table, the Fortran and the PL/M results 


are summarized. 


The PL/M results are the real internal values, not the 
displayed values. The Display routine gives only approximate 


values. 


This table shows for each run, the computed values of 
time of weapon release and corresponding aircraft's distance 
Wometarger and down range travel. It also shows the actuai 
values of distance to target and down range travel at the 


computed release time. 
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FORTRAN AND PL/i RESULTS: 


=n wee oe Se = —— oe ee ee eee 


Distance to target 
Down range travel 
fabulated values: 
Time 


Distance to target 


COHLUTED DORN RANGE TRAVEL 


Run #2 


= ee oe a= 


ee ques ewe 8 ee oe 


Tine 


Distance to target 
Dorn tange travel 
Tabulated values: 


ae wee ee Qe ame ame ang ae ae oo ome eS cS me ae 


Distance to target 


COMFUTED DOWN RANGE TRAVEL 


— yy oa —~ am 


Release point: 
Tine 

Distance to target 
Down range travel 
Fabulated values: 
Time 


Distance to target 


COMPUTED DOWN RANGE TRAVEL 


— wae oe —— = 


Distance to target 


Down Tange travel 


om eee eee eS 2 eee ee oe o_o oe ee SS oe 


Distance to target 


COMPUTED DOWN RANGE TRAVEL 


US 


POEL rai 


10.84 
ie dere 
1831.4 


10.84 
163125 
N85255 


11.14 
ee SieZ 
eo oe 


17.174 
1894.7 
eS. Weis 


11.18 
fea 7.3 
183743 


11.78 
1836270 
vos 1.2 


Vic 26 
2456.6 
2456.6 


pale 0 
DASE S10) 
2455.9 


PL/M 


1057 
1865.6 
W370. 


1027 
1864.3 
1888.0 


Lie © 
ell (a8 
1820.6 


Wales 
1337. 5 
(ce ore, 


11.0 
186456 
1861.3 


17.0 
1366.0 
1861.0 


i123 
2444.2 
2434.6 


Pies 
2HU46.14 
2466.4 





The results show two kinds of disagreements: between 
FORTRAN and PL/M results and also between the extrapolated 


and computed values of down range travel. 
Both disagreements are due to: 


1. The PL/M program that runs in the Ballistics conputer 
to calculate the down range travel of the weapons is 
Peeaveing results that differ from its FORTRAN counterpart. 
This fact was not observed before. It so happens that with 
the data used in this implementation the PL/M program was 


not aS accurate as it had been in previous tests. 


2. The extrapolation of down range travel useS computed 
values, which in spite of being originated from smoothed 
Gace CcOntain roundofrc errors. when a tied degree 
polynomial is fit to four of these points, the noise present 
disturbs the extrapolation. A second degree polynomiai 
through the same four points would probably give better 


results. 


The time the Executive computer took to find the release 
point was on the order of 0.9 seconds, which is somewhat 
iLicieileys This time can be reduced by using a faster 


extrapolation technique. 
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1. The hardware system used to demonstrate the tactical 


system is complete and working properly. 


2. A different extrapolation technigue is necessary in 
order to produce faster and more accurate results for the 


down range travel. 


Seine bdilestics computer is time critical in this 
system and pecause of that the Executive computer should not 


hold up the Ballistics ccmputer. 


This thesiS waS not able to prove that the implemented 
Systen would meet the requirements of a real tactical 
system. It merely provides a framework with which further 
tests could be exercised. These tests should include 
different methods of extrapolations and ali types of 
weapons. In particular, with respect to the extrapolation 
technigue, a straight line fit through the last three data 
points seemS to be promising both in accuracy and speed of 


computation. 
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APPENDIX A 


THE CONVERSION OF DATA PROGRAM AND OUTPUT 


This appendix shows the program that performed the 
conversion from the encoded values of range, elevation and 
Mermgtn Of alrcraft as tracked by radar into decimal 
representation used by the program to test the extrapolation 
routine and also to produce the data for the PL/M program 
which uses a three byte floating point number 


representation. 


Following the program iisting, the outputs are shown. 


The floating point number used is composed of three 
bytes, as follows: The two higher order bytes constitute 
the mantissa. The first bit of the mantissa is always. one, 
except in the case of the number zero, l1.e. the number is a 


normalized binary number. 


The lower order byte represents the exponent in tne 
following manner: The first bit is the mantissa Sign bit. 
If the number is negative, this bit 1s one. The following 
bits represent the exponent of two. This exponent is biased 
by hexadecimal forty, l.e. from this number and up, the 
exponent is positive and from this number down, the exponent 


is negative. 
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To make it clearer, some examples are given: 


PL/M floating point number decinal value 
E100 49 450.0 
AQ 00 C4 =55..0 
80 00 3F 0.25 
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APPENDIX B 


DoreeGALMArOLATING TESTING PROGRAM AND OUTPUT 


This appendix shows the program that was used to test 
the extrapclation routine on the System /360, together with 


mes Output. 


The data used in this program is the same used in the 


PL/M implementation. 
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/% 


CD9K,042H, : 

OB4H,O02FH,044H,0D8H,054H; 0C5H,OADH;, O08H,045H, OBEH, 

055th; OC5H, 

aC) ae 20E TH,OFTH,0C4H,005H,03AH,02DH, OEBH, 
c& i cre 

OABH, OAAH,03EH,O0C6H, OB1LH,OBEH,OCDH, O8EH,03FH,OA3H 

O90H,OBEH, ai 

O96H, 028H,03DH, OA6H,085H, OBSH; OAOH;,05AH,03 38H); 

DECLARE THE MACH CUT MATRIX Oa laerar d 

ARE CIVALUE (18) BYTE 

ITT AL{OD5H, O81H,040H; OFAH, O1CH ,040H »09FH, 0 32BH,040H; 0£2H, 
O8BFH,040H, 

084 ,018H,041H, OA 6H, 066H,041H) 3 
DECLARE THE WPNCODE MATRIX CONTAINING THE VARIABLE FOR 
EACH WEAPOK */ 

ARE WPNCOCE (585) BYTE 

WEAPCN CONSTANTS FOR THE MK 43 UNRETARDED */ 

ITLAL(4,000H, OOOH? OOOH; OA 7H, 02 7H} 038H; OAQHs000H, 04 3H, 
OCOH,000H,042H,; 

hHEAPCN CUCNASTANTS FOR THE MK 57 UNRETARDED */ 

45 000H,000H;, O00H, OCEH, 06 BH,0329H;0A0H ,900H,043H, 
OCOH,900H,042H, 

WEAPCN CONSTANTS FOR THE MK 61 UNRETARDED *«/ 
4,000H, 900i; OCOHs 083H;,066Hq 039H, OAOH ,000H19043H, 
OCOH,000H;0424h, 

WEAFCN CONSTANTS FOR THE MK 116 WETEYE */ 
21080H, O90H;: 039H, 084H,070H, 038H,CTOH,SOCH, 042H, 
O8CH, VOOHs, 0424, 

WEAFOCN CORSTANTS FOR THE MK 76 WITH LUG */ 
21080H,90CH ,039H,OD0H,090H;0325H,0COH ,00CH,042H, 
OCCH,900H,042H, 

WEAFOCN CONSTANTS FOR THE MK 77 FIRESOMB */ 
4¢9000H ,QOO0OH, OCOH, OAEH, 026H,038H;C380H,O000H;042H;, 
O80h, OOOH O41H, 

WEAPEN CONSTANTS FOR THE MK 81 ¥/ 
1,0A4H,O081Hs» 942H, COOH, GOIH ,OUGH;0AJQH,00CH;,643H, 
OC CH, OOOH, 042H, 

WEAPCN CONSTANTS FOR THE MK 81 SNAKEYE UNRETARDED *«/ 
ers oa OUM, OCOH,sO0A0H., 0 CoH, 034H,0CCH, O00CH, 042H; 
O8CH;, 000H,042H, 

HWeaeieN GUNSTANTS FOR THE MK 82 MECH FUZE */ 
1,084H,010H,042H, OO0H,0Q000H, OG0OrF;s0A0H +} 000H;043H, 
OCOH, OOOH, 042h; 

WEAP@N CUCASTANTS FOR THE AK 82 ELECT FUZE */ 
1,0BFH,O213H; O41H, OOOH, O00H,0COH,OA0OH ,OO00H,043H, 
OCCH,;0900H, 042k, 

Weriee@nmeewoianio FOR THE 4K 83 MECR FUZE */ 
1,0ABH,0EAH;: 0O41H, OOOH,000H, O000H; OAOH , OOOH, 04 3H, 
OCCH,000H, 0424; : 

WEAPON CUNSTANTS FOR THE MK 82 ELEC FUZE */ 
1,O9AH,CELH, O41LH, OOSH,CO0H,000H,0A0H;,000H,043H, 
OCOH,000H,042Hs 

WEAPON CONSTANTS FOR THE MK 84 */ 

LyO890H;000H; O41H, O00Hs OOOH »yOO00H, OAQH ,000H,043H, 
OCOH,000H; 042H, 

WEAFCN COASTANTS FOR THE MK LI? AL ¥/ 

1,0C 7H, OAEH;042H, OAOH, 04DH , 0B 7H,OAOH sO00H,C43H; 
_ OCCH;,;000H,04?H, 

WEAFCN CONSTANTS FOR THE MK 86 WET SAND FILLEC */ 
1,;00DFH,9D02H;,042H: 0O00H:000H;,000H,0COH,00CH;,042H, 
C8CH, O00H,042R, 

WEAPGN CONSTANTS FOR THE MK 88 WET SAND FILLED */ 
1,0CDH,0 70H, 041H; OOOH, GOOH, OOGH,0ADH -O900H,042H;, 
OC CK; 000H?7042H y 

WEAPGN CONSTANTS FOR THE MK 82 SNAKEYc UNRETARDED */ 
4,000H,DI0H, OOOH, 9FO0H,028H1039H:0COHs000H: 0O42Hy 
O8CH,O00H, 041H, 

WEAPCN CONSTANTS FOR THE MK 82 SNAKEYE RETARDED */ 
MeOOOH + OUCH, 000m, Or UH. 0 20h, 035H, OOOH ,OO00H,Q00H, 
O8CK,000H;,042Hy; 

13:2,000H;00GH, O00H; OSAH, 06 7H,03BH,0C2H,O8FH,03FH, 
OAFH, CCiH, O3S3EH. 
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